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Abstract 

We introduce a new parametrization for the dark energy, led by the same idea to the linear expansion of the equation 
of state in scale factor a and in redshift z, which diverges neither in the past nor future and contains the same number 
of degrees of freedom with the former two. We present constraints of the cosmological parameters using the most 
updated baryon acoustic oscillation (BAO) measurements along with cosmic microwave background (CMB) data and 
a recent reanalysis of Type la supernova (SN) data. This new parametrization allowed us to carry out successive 
observational analyses by decreasing its degrees of freedom systematically until ending up with a dynamical dark 
energy model that has the same number of parameters with ACDM. We found that the dark energy source with a dy¬ 
namical equation of state parameter equal —2/3 at the early universe and —1 today fits the data slightly better than A. 

Keywords: dark energy experiments • dark energy theory [] 

1 Introduction 

The recent high precision data is in very good agreement with the six parameter base ACDM cosmology mm- However, 
we still cannot consider it as the final cosmological model. It may represent, rather, a very good approximation or a 
limiting case of a more general theory, which does not necessarily reduce to vacuum energy in terms of a cosmological 
constant A and may deviate from A considerably in the far future and/or past. Two well known problems related with 
the A assumption, the so called coincidence and fine-tuning problems, may be signs for a dynamical nature of dark 
energy [3l |4l O |6l [3 E] . Some tensions between the A assumption and high precision data could be resolved in case of 
evolving dark energy (DE) (see for instance |9l llOi fTT] and references therein). 

Using the one parameter extension to the ACDM, wCDM adopts a spatially flat universe and constant equation of 
state (EoS) parameter w for DE, Planck collaboration [5] gives w = —I.ISIq'js from combined Planck+WP-|-highL+BAO 
data and BOSS collaboration uni gives w = —0.97 ± 0.05 from the most recent combined Planck-|-BAO-|-CMB data. 
Constraining the possible evolution of w on the other hand is difficult even with the powerful BAO, SN and CMB data 
that range at different redshift values. Hence, probing the possible DE evolution usually involves the introduction of 
a phenomenological parametrization for its EoS parameter involving a couple of free parameters. The most widely 
used parametrization for constraining evolving DE is the linear expansion in scale factor a which is known as CPL 
(Chevallier-Polarski-Linder |12[ 113) 1 parametrization: w{a) = too + Wa{l — a), where wo and Wa are real constants. 

It has been first proposed as an alternative to the linear redshift parametrization of DE: w(z) = wo -1- WzZ, which in 
contrast to CPL grows increasingly unsuitable at redshifts 2 3> 1 and hence cannot be used with high-redshift data, 
e.g., CMB. However, even the tightest constraints on the parameters of CPL from combined Planck-l-BAO-l-SN allow 
a quite flexible range for the evolution of the EoS parameter for DE m- Besides this, compared to ACDM, the 
improvement in the success of fit to the data using this parametrization is not significant [3 HOI lll | . These might be 
signaling for that CPL is not an adequate choice for describing the possible evolution of dark energy. It is already 
mentioned in many studies that the divergent behavior of CPL not only prevents one to make plausible predictions on 
the future of the universe but also shows that it cannot genuinely cover theoretical models of dark energy. Throughout 
the literature there have been various DE parametrizations introduced by considering different reasoning and strategies 
for obtaining mathematically well behaving and physically acceptable models to obtain tighter constraints on DE and 
make plausible predictions on the future of the universe [nilT^fT51fT71fT51[T^[^[^l^l^[M[^[^[?7| - 

In this paper we shall use a new parametrization for describing the DE source that may be taken as a natural 
extension of the same idea that give rise to linear redshift and CPL parametrizations that correspond to the first two 
terms of Taylor expansion in z and a, respectively. We start by introducing a fluid with an EoS parameter linear in 
time t. Then, we write its EoS parameter in terms of the scale factor a by considering general relativity and utilize the 
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obtained w{a) for describing the DE component of the universe in our observational analyses. We carry out successive 
observational analyses by decreasing the number of degrees of freedom (DoF) until we end up with a dynamical DE 
model having no additional parameters compared to ACDM. 


2 Dark energy parametrization 


In this section we first derive an EoS parameter in terms of scale factor a, which yields linear EoS parameter in time 
t when the universe is filled only with this fluid (here tilde denotes a universe contains only the fluid described by 
the EoS linear in time t). We then in Sec. 2.2 make use of this form of w{a) for describing the DE source in the 
physical universe that contains sources other than DE also, by adopting d —>■ a, where a represents the scale factor of 
the physical universe. 


2.1 Equation of state parameter linear in time 

Let us consider a fluid described by an EoS parameter expressed as a first order Taylor expansion in time t, that is: 


w = Wo + wi(l - i), 


( 2 . 1 ) 


where Wq and Wi are real constants and t > 0 is the normalized time. It should be noted here that w would not 
diverge provided that the time is restricted as in the big ri]:Q cosmologies. We construct a relation between the fluid 
described by the EoS parameter linear in time (2.11 and the scale factor by considering the general theory of relativity. 


Accordingly, we consider the Friedmann equations for spatially flat RW spacetime in the presence of a single fluid 
-2 ~2 ~ 

written as follows: 3^ = Kp and Sj + 2| = —Kp = —npw, where a dot denotes derivative with respect to time t. 
Eliminating p between these two equations and then using (2.11 we get the corresponding deceleration parametei]^ 


dd 1 Up /i 

1 = = 2 2 Fo+w'uI-dJ ’ 


whose solution, for wi ^ 0, reads 


a = ai exp 


arctanh 


piJlt—(- u;q+p<;i + 1) 

\/ {mo+wi+^^—ciwi 


\J (wo + Wl + 1)2 — ClWl 


( 2 . 2 ) 


(2.3) 


where 5i and ci are integration constants. For the case Wi — 0, we would get the standard power law solution as 
5 oc p/3(i+“o) ^ith p oc We demand a big bang a = 0 at t = 0, and expanding universe a > 0 for t > 0 and 

5 —>■ +00 at a finite time, say, as t —>■ Ibr where Ibr is the big rip time, so the EoS parameter never diverges. One may 
check that the first condition is satisfied by choosing Ci = 0, and in addition to this the latter two conditions imply 


the following relations 2 + 2wo + 2wi > Wi > 0. Accordingly, substituting ci = 0 in (2.31, our solution reduces to 


d = di exp 


/I arctanh . , , , 

3 Wo + wi + 1 


- 1 


with 2 + 2wo + 2wi > wi > 0. Isolating t in this solution by setting d{t = 1) = 1 we obtain 

2(1 + Wo + wi) 


t{d) = 


wi + (2 + 2wo + wi)a 2 


Finally using tid) in (2.11 we obtain the following EoS in scale factor d 

2(1 + Wo + Wi) 


w(a) = Wo + wi 


1 - ^ 


Wi + (2 + 2wo + wi)c 


-3(l + mn + u,i) 


(2.4) 


(2.5) 


( 2 . 6 ) 


*See 1281 for big rip cosmology. 

tit may be noteworthy that a deceleration parameter linear in time for the observed universe was obtained from a higher 
dimensional cosmological model in dilaton gravity in [29], where its value at f = 0 is depending on the number of extra dimensions. 
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with 2 + 2wo + 2wi > Wi > 0. Next, using the energy-momentum conservation equation we obtain the energy density 

P- 




1 + Wl 


(™0+™i+l) _ 1 
‘2wo 2mi “t" 2 


(2.7) 


where po = p(2 = 1). We notice that 


—3(iyo + l) ln(a)—4i[;o 
3(tf;o + l) ln(a)—4 


and p —>• po — 1 -I- 


3(mo + l) 


■ ln(a)j as ici —>■ — 1 — wo- 

The EoS parameter (2.61 is bounded both as a —>■ 0 and a —>■ oo, namely, w —>■ too + as a —>■ 0 and w —>■ 
— Wo — ici — 1 as a —>■ oo. We note that a —>■ 0 as t —>■ 0 and 2 —>■ oo as t —>■ 2(1 -|- Wo + wi)/wi, hence the time could 
not take arbitrarily large values unless Wi = 0, which also explains why EoS parameter linear in t does not diverge in 
future in this case. 

The recent high precision data is in very good agreement with A as the DE source in the universe. However, it 
may represent, rather, a very good approximation or a limiting case of a more general theory, which does not reduce 
necessarily to vacuum energy in terms of a cosmological constant. In this regard, being Wde ~ —1 is in very good 
agreement with precision data, it seems physically more acceptable to consider an approximate description of possibly 
dynamical DE with a parametrization that is mathematically well behaved and may never depart from Wde ~ — 1 a 
lot. We note that the EoS parameter we derived has these properties and hence may be adopted for parametrizing the 
EoS of the DE source in the universe containing other sources also. 


2.2 Parametrization of dark energy source 

The EoS derived above was hrst given in m starting from the linear deceleration parameter in time using general 
theory of relativity and then utilized in [SlISl] for describing the effective EoS parameter averaging all the ingredients 
of the universe to constrain particularly the kinematical properties of the expansion of the universe. In this work, on 
the other hand, we follow a quite different method and adopt the EoS written in terms of 2 for describing only the 


DE component of the universe, rather than the average ingredient of the universe. The cost of this is that (2.6 1 will 


not correspond to the linear EoS in physical time t in the actual universe anymore, unless the universe is dominated 
by this fluid. The reason being is that including the other sources such as cold dark matter (CDM) in addition to DE 
the evolution of the scale factor will deviate from the one given in (2.41, i.e., the relation between the scale factor a 
and time t will be different, until the DE source becomes dominant over all other sources in the universe. 

Let us now adopt the EoS parameter given in (2.6 1 for describing DE source by setting 2 —>■ a and then write it in 
a more useful form by defining Wi = wo -|- wi: 


Wde = Wi - ' 


2(wi — wo)(l -I- Wi) 


(wi — Wo) -I- (2 -b Wi -I- Wo) ( 


which also recasts the energy density (2.71 as 


(0) 

pde — Pde 


Wi - Wo 3 (. 

-a*' 

2wi -1-2 


„i+i) Wi -I- wq -b 2 3 

2wi -b 2 


i+i) 


( 2 . 8 ) 


(2.9) 


where a is the physical scale factor, and the conditions on the parameters are now given by —2 — Wi < wo < Wi. Our 
parametrization has the same number of degrees of freedom with the CPL but while the variation of the EoS of DE 
with respect the scale factor is constant in scale factor w^pl = —Wa in CPL, here it is a dynamical quantity: 


^de 


3(wi -b l)"[(wo -b If - (wi -b 

|^(wo — Wi)ai“’‘+3 — (wo -b Wi -b 2)j 


( 2 . 10 ) 


where a prime denotes derivative with respect to the scale factor a. 


We note that the form of the energy density (2.91 we obtained by the reasoning we followed in the previous section 


has a symmetry for Wi -b 1 — >■ — Wi — 1 and hence it is enough if we consider the cases Wi > — 1 only. In accordance with 
this, we see from (2.8 1 that Wde — >■ Wi as a —> 0 (a —>■ oo) and hence Wi denotes the initial value of the EoS parameter 
of DE. Similarly Wde —t — Wi — 2 as a —^ oo (z —>■ — 1), and hence we can write the final value of the EoS parameter as 
Wf = —Wi — 2. We also note that w'd^ can take only negative values as long as we stick to the range —Wi — 2 < wo < Wi 
for the parameters, but take positive values if this condition could be violated. The radiation and pressure-less matter 

4 2 

components of the universe enforce universe to expand as a oc 12 and a oc 13 , respectively. The DE source described 


by (2.91, on the other hand, enforces universe to expand according to the kinematics we discussed in the previous 
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section. In the physical universe, say, in the presence of all these components, the expansion of the nniverse will be 
much complicated as it will be determined by the joint effect of these components. We also note that the conditions 
on the parameters of the hypothetical fluid, we are now utilizing for describing DE, were initially considered to obtain 
a physically reasonable expansion history for the universe filled by this fluid only. However the presence of radiation 
and pressure-less matter in addition to a DE source described by (2.91 can lead to viable expansion histories for the 
universe even if the parameters Wi and wq take values out of the range —Wi — 2 < wq < Wi and hence provides us 
freedom to violate this condition. Accordingly, in the next section we will carry out observational analyses without 
taking this condition into account, so that our analyses will not be restricted to the cases ui^e < 0. However we will see 
below that the concluding values for Wi and wo upon the three successive observational analyses, where we decrease 
the number of degrees of freedom systematically, satisfy the condition —Wi — 2 < wo < Wi. 


3 Observational constraints 


We consider a spatially flat RW spacetime and the Hubble parameter H in the presence of radiation Hr, pressure-less 
fluid/(CDM-|-matter) Hm and the dark energy Hde described by (2.91 is given by 

Wi -I- Wo -I- 2 




- o(0) 




2wi -b 2 ^ ^ 2wi -b 2 


(1 + ^) 




(3.1) 


where the density parameter of the DE = 1 — and Wi > —1. Here subscript/superscript 0 indicates 

the values of the parameters today. The density parameter of radiation is = 2.469 x 10~^h~^(l + 0.2271Vefr), 
where h = ii/o/100kms~^Mpc“^ [I]. We consider a model with standard matter and radiation content, including three 
neutrino species with minimum allowed mass = 0.06 eV. Throughout the analysis we assume flat priors over 

our sampling parameters: wo = [—2.0, —0.5], wi = [—1.0, 2.0] for the dark energy EoS parameters, = ]0.05,1] for 
the pressure-less matter density parameter today, = [0.02,0.025] for the baryon density today and h = [0.4,1.0] 

for the reduced Hubble constant. It might be noteworthy to note that the radiation energy density Hq 

today is not subject to our analysis since it is well constrained, such that it has a simple pr = ^Tqmb relation 
with the CMB monopole temperature (see [33] for further details), which today is very precisely measured to be 
^CMB = 2-7255 ± 0.0006 K [31]. 

In order to perform the parameter space exploration, we make use of a modihed version of a simple and fast Markov 
Chain Monte Carlo (MCMC) code that computes expansion rates and distances from the Friedmann equation, named 
SimpleMC. This code already contains a compressed version of the Planck data, a recent reanalysis of Type la supernova 
(SN) data, and high-precision BAO measurements at different redshifts up to 2 < 2.36 |10| . 

In three successive subsections we carry out observational analyses: i) wo and Wi are both free parameters: DoF 
are the same with CPL, ii) wo is free but Wi is fixed to certain values, so that DoF are the same with wCDM, and iii) 
both Wo and Wi are fixed to certain values so that DoF are the same with ACDM. 


3.1 General case: free wq and w\ 

We first constrain the cosmological parameters by taking both wo and Wi as free parameters so that we have the same 
number of DoF with the CPL parametrization. We summarize the results for the data sets Planck-bBAO, Planck-bSN 
and Planck+BAO-bSN in Table In Figure we show 1-D and 2-D posterior distributions of wq and Wi for the 
following data sets Planck-bBAO, Planck-bSN and Planck-bBAO-bSN. In Figurej^we give 2-D posterior distributions 
for Wo — h and wo — . 

We observe that the matter density parameter and the reduced Hubble constant h are obtained almost the same 
in three combinations of data sets and Planck+BAO-bSN data set leads to the tightest constraints; = 0.304 + 0.009 
and h = 0.677+ 0.011. 

In all cases the best-fit parameters are consistent with ACDM model, i.e., (wo,Wi) = (-1,-1). We observe that 
Planck+SN data lead to a tighter constraint on wq while Planck+BAO data lead a tighter constraint on Wi. The 
central value of wo is almost equal to —1 from both Planck+BAO, Planck+SN and Planck+BAO+SN and its value 
from the different data combinations changes by less than la. Using three data sets at the same time we obtain 
Wo = —0.99 + 0.06 (Planck+BAO+SN) with la having the tightest constraint on wo with a central value almost 
exactly equal to —1. Using CPL, on the other hand, wo is not constrained very well and the central value of wo 
is obtained not only quite higher than —1 but also significantly different from three different data combinations; 
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Table 1: Parameter constraints. SN data cannot constraint Wi as good as wq while it is the other way around for BAO data. 
For two-tailed distributions the results are given in Icr and for one-tailed distributions given in 2c7. 


Data set 


h 

Wo 

— 1 < Wi 

Planck±BAO 

0.306+g:°i^ 

D f;7Q+0-020 

U.D/ O_o.oi8 

-0.97 ±0.10 

< -0.30 

Planck±SN 

U.OU^_0.oi9 

0.677 ±0.019 

_ 02+0.07 

-L-UZ_o 08 

< -0.27 

Planck±BAO±SN 

0.304 ± 0.009 

0.677 ±0.011 

-0.99 ±0.06 

< -0.42 



Figure 1: 1-D and 2-D posterior distributions for the parameters wq and W[; data combinations are indicated in the legends. 




Figure 2: 


2-D posterior distributions for wq — h and H 


( 0 ). 
m ) 


data combinations are indicated in the legends. 


Wo = —0.58 ± 0.24 (Planck+BAO), wo = —0.90 ±0.16 (Planck±SN) and wo = —0.93 ±0.11 (Planck±BAO±SN) |10j . 
Constraint on Wi from Planck±SN is quite loose such that — 1 < tCi < —0.27 (2(j); so that the DE could either behave 
like cosmic strings or a cosmological constant at very high red-shift values. Planck±BAO data constrain Wi better 
and constrict the allowed range to — 1 < tCi < —0.30 (2a). The tightest constraint is obtained from Planck±BAO±SN 
data; —1 < Wi < —0.42 (2ct). According to this initially dark matter and cosmic strings like DE scenarios are not 
viable. Considering the constraints on wq and Wi together we observe that A as the DE source is doing great but 
even considerably large deviations from A are still allowed, for instance, (wo,Wi) = (—1, —|) is also perfectly allowed 
in this picture. The best fitting model to Planck±BAO±SN data has = 46.63, representing an improvement of 
Ax^ = 0.37 compared to ACDM for which Xacdm = 47.00, for the addition of two extra parameters. This is an 
improvement almost the same with the one in the CPL model Axcpl = 0-42 uni. Accordingly in terms of information 
criteria the improvement is not sufficient enough to justify the addition of two extra degrees of freedom either in our 
parametrization or CPL parametrization and additionally there is no reason to prefer one over the other among these 
two parametrizations. 
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3.2 Free wq and fixed Wi 

We may try to get more information nsing our parametrization by redncing its DoF to the that of the wCDM 
parametrization by fixing Wi. We now carry out the observational analyses nsing Planck+BAO+SN data by setting 
either w; = 0 (DE starts like pressure-less matter), Wi = — | (DE starts like cosmic strings) or Wi = — | (DE starts like 
cosmic domain walls). In Tablej^ we snmmarize the results, including the minimum values, from Planck-|-BAO+SN 
data. For a comparison, the constraint on wCDM parametrization from Planck-|-BAO-|-SN data is wq = —0.97 ± 0.08 
(Icr) [TO]. We give the 1-D probability distributions of wo in our model in Figure]^ We note that the central value 
of Wo shifts from the values less than —1 to the values higher than —1 and the x^ improves as Wi goes from 0 to — 1. 
Given that the case Wi = 0 yields very large x^ relative to other three, it can be rejected. The most interesting point 
in Table is that wo takes the closest value to —1 in the case Wi = — | rather than in the case Wi = — 1 and these two 
cases have almost the same x^ values that are considerably lower than the cases Wi = 0 and Wi = — This may be 
interpreted as a signal for a DE with a dynamical EoS parameter starting from values ~ — | at earlier times and then 
approaching a cosmological constant as the universe expands, as noted also by m- 

Table 2: Mean values with la errors of the parameters of the model for the particular values = 0, Wi = = —^ and 

Wi = —1 from Planck+BAO+SN data. 



h 

Wo 

Wi 


0.289 ±0.008 

0.688 ±0.010 

1 1 on+0.028 
l.ibU_Q Q25 

0 

54.04 

0.298 ±0.008 

0.680 ±0.010 

i.UDO_o 035 

1 

3 

48.44 

0.303 ±0.009 

0.677+°?? 

U.yi^O_o 042 

2 

3 

46.90 

0.305 ± 0.009 

0.676+°l? 

-0.966 ±0.053 

-1 

46.64 


1.0 

0.8 
I 0.6 

0.4 

0.2 

0.0 


Figure 3: 1-D posterior distributions for wq for the particular values tcj = 0 (blue), Wi = —^ (green), Wi = —| (red) and 
lUi = — 1 (black) from Planck+BAO+SN data. 


3.3 Fixed wq and fixed Wi 

We observe from Table that wo persists on yielding values around —1 in spite of the large differences between the 
fixed values of Wi. Hence we now fix also wo = —1 along with Wi = 0, Wi = —Wi = —| or Wi = —1. Doing 
so our model now yields the same number of DoF with the ACDM model and the latter case {wQ,Wi) = (-1,-1) 
corresponds to A while the former three cases correspond to dynamical DE models. We summarize the constraints 
from Planck-l-BAO+SN data set with their minimum x^ values in Table[^ We note that the case {wo,Wi) = (—1,0) 
and {wo,Wi) = (—1, —|) should be ruled out due to their very large x^ values. We note on the other hand that the 
cases (wo, wi) = (—1, — |) and (wo, wi) = (-1,-1), i.e. ACDM, have almost the same low x^ values; X^i -2/3) = 46.88 
and Xacdm = 47.00, respectively. However, it is noteworthy to note the case {wo,Wi) = (—1,—|) has the lowest x^ 
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Table 3: Mean values with Icr errors of the parameters of the model with wq = —1 for the particular values Wi = 0 . Wi = — ^ 

9 0 

and nji = — g and values from Planck+BAO+SN data. The case = —1 corresponds to ACDM. 


4 ijji 

h 

Wo 

Wi 


U.diO_o.o05 

0.605 ± 0.002 

-1 

0 

173.34 

U.OUJ_Q Q08 

0.664 ± 0.005 

-1 

1 

3 

51.80 

U.OUZ_Q QQg 

0.678 ± 0.006 

-1 

2 

3 

46.88 

0.302 ±0.008 

0.682 ± 0.006 

-1 

-1 

47.00 


value with a difference over the ACDM given by Xacdm ~ X^-i - 2 / 3 ) = 0.12 , although not statistically significant. 
In this case we find from (2.101 that < 0 through the history of universe and w^g(z ~ 0) = today, while 
w' is always null in the case of A. According to this we have a dynamical DE model that fits to data equally better 
than A but yet does never cover/mimic A. We plot Wde and dwde/dz = —Wde/(1 + versus redshift 2 in Fig. [^for 
demonstrating how the dynamics of DE in case {wo,Wi) = (—1, — |) deviate from A. 

Two key parameters in cosmology are the deceleration and jerk parameters that are defined as q = — = 

^ ^ J ~ afi^ = 9 + 2(7^ + (1 + 2 ) ^, respectively. The negative valnes of deceleration parameter imply that 

the universe is expanding with an accelerating rate and values less than —1 indicate a super-exponential expansion. 
The jerk parameter on the other hand is a very useful parameter for investigating the deviation of a cosmological model 
from ACDM since its value stays pegged to unity in ACDM (ignoring the contributions other than A and pressure-less 
matter), while it is in general a dynamical quantity in onr model. In the case {wo,Wi) = (—1,—|) the deceleration 
parameter evolves from 0.5 as in the ACDM and goes ever monotonically to — | which indicates a big rip end of the 
universe while the universe approaches asymptotically de Sitter universe with a deceleration parameter equal to —1 in 
ACDM model. We depict the evolution of q and j in redshift 2 for (wq, Wi) = (—1, —|) and ACDM in Figj^using the 
values with Icr errors from Table|^ We find that the current value of the deceleration parameter is qo = —0.545±0.014 
(1(t) and the universe starts its accelerating expansion at redshift value Ztr = O.OdTio'Qis (1'^) while these values are 
qo = —0.547 ± 0.012 (Icr) and 2tr = 0.665 ± 0.021 (Icr) in ACDM model. We observe that the jerk parameter exhibits 
a non-monotonic behavior but stays at values very close to unity, and find its current value is j{z = 0) = 1.087 ±0.001 
(Icr). Although the expansion history of the universe is almost the same in both models, they predict completely 
different futures. In the case (wo,Wi) = ( —1,—|) the universe enters into a super-acceleration regime (q < —1) at 
redshift 2s = ~0.487jjQ'gQ| (Ict) and ends with a big rip. 




Figure 4: and d?rde/d 2 in terms of redshift 2 for the case (wo,Wi) = (— 1,— |) (solid lines). The dotted lines correspond 

to the case cosmological constant. 


4 Concluding remarks 

The new parametrization of DE we introduced in this paper allowed ns to carry ont snccessive observational analyses 
by decreasing its degrees of freedom systematically until ending up with a dynamical dark energy model having no 
additional parameters compared to ACDM. We first fixed the initial value of the EoS to some particular values and 
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Figure 5: Deceleration parameter q and jerk parameter j versus redshift z for the case (wo,Wi) = (—1, —|) (solid lines) and 
ACDM (dotted lines). 


found some indications in favor of an evolving DE. We then additionally fixed the EoS parameter to —1 for the present 
epoch of the universe and found that EoS parameter that equals — § at high redshifts and —1 at the present epoch 
gives a slightly better fit to the combined Planck+BAO+SN data than A. The variation of the EoS parameter in 
scale factor, on the other hand, is dynamical but always less than zero, which means that even today, when its EoS 
parameter is equal to —1, the DE does not genuinely mimic A. Thus we gave a cosmological model that can compete 
with ACDM but involving a dynamical DE. We also note that DE in this model passes below the phantom divide line 
in the future, hence leading to a big rip end of the universe. Accordingly, because the universe lives for finite time and 
the EoS parameter of DE evolves around —1, hence the energy density of DE is also dynamical, throughout the history 
of universe, our model relieves also the problems related with cosmological constant assumption of ACDM model. 
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